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We consider oscillations of the low energy (sub-GeV sample) atmospheric neutrinos in the three 
neutrino context. We present the semi-analytic study of the neutrino evolution and calculate char- 
acteristics of the e-like events (total number, energy spectra and zenith angle distributions) in the 
presence of oscillations. At low energies there are three different contributions to the number of 
events: the LMA contribution (from ^-oscillations driven by the solar oscillation parameters), the 
!7 e 3-contribution proportional to S13, and the U e 3 - induced interference of the two amplitudes driven 
by the solar oscillation parameters. The interference term is sensitive to the CP- violation phase. We 
describe in details properties of these contributions. We find that the LMA, the interference and U e 3 
contributions can reach 5 - 6% , 2 - 3% and 1 - 2 % correspondingly. An existence of the significant 
(> 3 — 5%) excess of the e-like events in the sub-GeV sample and the absence of the excess in the 
multi-GeV range testifies for deviation of the 2-3 mixing from maximum. We consider a possibility 
to measure the deviation as well as the CP- violation phase in future atmospheric neutrino studies. 
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Cn ' I. INTRODUCTION 

o 

^ The existing results on atmospheric neutrinos Q are well described in terms of pure „ „ u T oscillations 

with maximal or close to maximal mixing. Analysis of the Super-Kamiokande data gives [HQ] mass squared 
Oh, difference and mixing in the interval: 

Amj 2 = (1.3 - 3.0) x l(T 3 eV 2 , sin 2 2(9 23 > 0.9, (90% C.L.). (1) 

The results of SOUDAN Q and MACRO Q experiments are in a good agreement with The oscillation 
/\ ' interpretation Q has been further confirmed by the results of the K2K experiment [I| . 

H ' 

Till now no compelling evidence of oscillations of the atmospheric v e has been obtained. The Zv global 
analysis of the atmospheric neutrino data is in agreement with no i/ e -oscillations |(|. The best fit point 
coincides with zero sin #13 within la 

At the same time, after the first KamLAND result Q| we can definitely say that the j/ e -oscillations 
of atmospheric neutrinos should appear at some level. Indeed, KamLAND has confirmed the large mix- 
ing MSW (LMA-MSW) solution of the solar neutrino problem. The combined analysis of the solar and 
KamLAND data leads to values of the oscillation parameters Q : 

Am 2 ,! = (5 - 10) • 10~ 5 eV 2 , tan 2 12 = 0.3 - 0.5 . (2) 

The parameters ©, which we will call the LMA parameters, should lead to oscillations of the v e component 
in the atmospheric neutrino flux. 
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The oscillations driven by the LMA parameters have been discussed before fiot fTll fT^. f]~3L fl^ . ll5L TT^ | - 

It was marked in Ref. that the effect of sub-leading oscillations driven by Am^ is significant only for 
the sub-GeV events and the size of effect is at the level of the statistical errors. In Ref. [13 it was argued 
that the excess of e-like events in the sub-GeV sample favors the large mixing MSW solution of the solar 
neutrino problem. 

The detailed study of the effect has been performed in our previous paper , where it was shown that 
the neutrino oscillations with LMA parameters can lead to an observable (up to 10 - 12 %) excess of the 
e-like events in the sub-GeV atmospheric neutrino sample. The excess has a weak zenith angle dependence 
in the low energy part of the sample and a strong zenith angle dependence in the high energy part. The 
excess rapidly decreases with energy of neutrinos, and it is strongly suppressed in the multi-GeV range. 
These signatures allow one to disentangle the effect of oscillations due to solar Am 2 from other possible 
explanations of the excess. 

It was shown that the relative excess is determined by the two neutrino transition probability Pi and 
the "screening" factor: 

|i-l = P 2 (r cos 2 #23-1) , (3) 

where F e and F® are the electron neutrino fluxes with and without oscillations, r is the ratio of the original 
muon and electron neutrino fluxes. The screening factor (in brackets) is related to existence of both the 
electron and muon neutrino in the original atmospheric neutrino flux. The appearance of excess (or defi- 
ciency) depends strongly on deviation of the — v T (or 2-3) mixing responsible for the dominant mode 
of the atmospheric neutrino oscillations from maximal value. Indeed, in the sub-GeV region r *=s 2, so that 
the screening factor is very small when the — v T mixing is maximal. Due to this factor the excess is in 
general small even though the 2v— probability can be of the order 1. The probability P2, and consequently 
the excess, increase rapidly with Am 2 !. 

As far as the experimental results are concerned, there is a hint that some excess of the e— like events 
indeed exists in the sub-GeV range. Furthermore, the excess increases with decrease of energy within the 
sample |l8j . In comparison with predictions based on the atmospheric neutrino flux from ref. |l7j the excess 
is about (12 - 15)% in the low energy part of the sub-GeV sample (p < 0.4 GeV, where p is the momentum 
of lepton). It has no significant zenith angle dependence. In the higher energy part of the sub-GeV sample 
(p > 0.4 GeV) the excess is about 5% , and practically there is no excess in the multi-GeV region (p > 1.33 
GeV). 

The excess is within estimated 20% uncertainty in the original atmospheric neutrino flux. The analysis 
of data with free overall normalization leads to the best fit e— like si gna l which practically excludes the 
excess 0, 0, 0|. However, the recent data on primary cosmic rays |2Ct l2l| as well as new calculations 
of the atmospheric neutrino fluxes 0| change the situation. New results imply lower neutrino flux, and 
therefore larger excess which is difficult to explain by change of normalization |l8l Il9| . 

The v e oscillations can be also induced by non-zero 1-3 mixing and Ami responsible for the dominant 
mode of the atmospheric neutrino oscillations 0, E |H ||l |H |2fl |U El ill OH EH]. These oscillations 
require non-zero value of mixing matrix element U e $ . They are reduced to the vacuum oscillations for the 
sub-GeV sample. For the multi-GeV sample the Earth matter effect becomes important which can enhance 
the oscillations For neutrinos which cross the core, the dominant effect is the parametric enhancement 
of oscillations. The size of effect is restricted by the CHOOZ bound on U e 3 |32| . 

In this paper we will further study the oscillation effects driven by the LMA oscillation parameters using 
an updated experimental information. We will consider an additional effect in the sub-GeV sample induced 
by non-zero 1-3 mixing. We study effects of the interplay of oscillations with the LMA parameters and 
non-zero U e ^. In particular, we will discuss the interference induced by non-zero sin $13. Some preliminary 
results of this study have been published in |33| . The interference term depends on the CP-violation phase. 
We calculate effects of CP-violating phase and estimate a possibility to observe it in future atmospheric 
neutrino experiments. 
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The paper is organized as follows. In Sec. II we present the semi-analytical study of evolution of 
2>v system at low neutrino energies which correspond to the sub-GeV sample of events. We present the 
2>v neutrino transition probabilities in terms of the two neutrino probabilities. We calculate the latter 
numerically and study their properties. In Sec. II. C we give general expression for the flux of electron 
neutrinos. In Sec. III. we calculate the number of e-like events in the water Cherenkov detectors with 
and without oscillations. We consider the zenith angle and energy distributions of these events. We study 
separately effects of the LMA-oscillations (Sec. III. A), the U e 3 induced interference (Sec. III.B) and the 
CP-violation (Sec. III.C). In Sec. IV we discuss a possibility to measure the deviation of 2 - 3 mixing from 
maximum as well as the CP-violating phase in future atmospheric neutrino experiments. Discussion of the 
results and conclusions are given in sect. V. 



II. EVOLUTION OF THE NEUTRINO SYSTEM 



We consider the three-flavor neutrino system with hierarchical mass squared differences: Am^i = 
Anig << Anijj = Am 2 tm (see Eqs. JTEJ). The evolution of the neutrino vector, Vf = (v e , v^, i> T ) T ', is 
described by the equation 

.du f [UM 2 U^ -\ 

=HB- +v r'' (4) 

where E is the neutrino energy, M 2 = diag(0, Am^, Am^) is the diagonal matrix of neutrino mass squared 
eigenvalues, V = diag(V,0,0) is the matrix of matter-induced neutrino potentials with V = ^/2GfN e , Gf 
and N e being the Fermi constant and the electron number density respectively. The mixing matrix U is 
defined through Vf — Uv mass , where v mass = {vi,v 2 ,v 3 ) T is the vector of neutrino mass eigenstates. The 
matrix can be parameterized as 

U = U 23 D 5cp U 13 U 12 , D Scp = diag(l, l,e i5 ° p ), (5) 

where Uij = U{j (fty ) performs the rotation in the ij- plane on the angle 6ij . This parameterization coincides 
with the standard one up to (unphysical) renormalization of the mass eigenstates: v rnass = D* s v' mass . 



A. Propagation basis 



The dynamics of oscillations is simplified in the "propagation" basis v — (24, z>2, V3) T , which is related 
to the flavor basis by the unitary transformation 

v f = UD, (6) 

and the matrix U can be introduced in the following way. First, let us perform the rotation Vj = 
UzsDsQpUxzv' . Using Eq. we fi n d that the instantaneous Hamiltonian for the new states v' takes 
the form 



or explicitly 



H' = ^U 12 M 2 U\ 2 + U1 3 VU 13 



s : f 2 Aml 1 /2E + Vc 2 3 s 12 c 12 Am 2 . 1 /2E 
H' = [ s 12 c 12 Am 2 1 /2E cf 2 Am| 1 /2E ) . (7) 

V e Si 3 Ci3 
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(ci2 = cos^i2, S12 = sin^i2, etc.). Since the transformation U2 3 D$ cp Ui 3 is constant (no dependence on 
density and therefore time), the evolution equation for v 1 is given by the Schrodinger equation with the 
Hamiltonian H' . Notice that the CP-violation phase is removed from the equation for v' , and it appears in 
the projection of the flavor basis on v' only. So, the evolution of system is CP-symmetric. 

Second, let us make an additional (y' e — v' T ) rotation, v' = U[ 3 C>, which removes the off-diagonal terms 
H' 13 , H' 31 of the Hamiltonian Q). The angle of this rotation depends on the density: 

1EV 

tan(A6»i 3 ) « s 13 c 13 — — — . (8) 
Am 13 

Since the low energy part of the sub-GeV sample is produced by neutrinos with energies (0.1 - 0.5) GeV 
the factor in Eq. JHJ) can be evaluated as 

^=5.i. 10 -^.-g-. i^y <o.i 

Amf 3 g/cc lGev Am z 

for the matter density in the Earth (3 — 10) g/cc. That is, the additional rotation is much smaller than the 
vacuum (1-3) rotation. The angle 6™ 3 can be considered as a small correction to #13: 



013 = 013 + A9 13 « 013 1 + -r— 5- . (9) 
V Am 31/ 

Basically, $13 is the mixing angle in matter in the two neutrino approach: 6*13 w 8" 3 , and plus sign in Eq. J5J 
reflects the fact that matter enhances the mixing in the neutrino channel for the normal mass hierarchy. 
In the antineutrino channel the sign is negative. For the inverted mass hierarchy the sign of Am| x should 
change, and the correction to S13 is positive in the antineutrino channel and negative in the neutrino 
channel. 

As a consequence of the additional rotation we find the following. 

1. The correction to H'n element appears which can be considered as the correction to the potential: 

AH n =AV^ ~V S j 3 (l - %?V\ . (10) 



Am 13 



It can be safely neglected. 

2. The H23 and H 32 elements are generated: 



H 23 = H32 = S12C12S13^ ( ) (H) 



'"'13, 

which are again negligible. Also corrections to the H 33 element can be neglected. 

Combining the rotations introduced above we find the projection matrix U © which defines the 
propagation basis: 

Cl3 S13 

U = U 23 Ds OP U 13 = U 23 Ds OP U 13 U{ 3 = \ -s 13 s 2 3e iScp c 23 5i 3 s 2 3e 1 ^ | . (12) 

-si3C 2 3e lScp S23 C 13 C23e lScp 

In the propagation basis the evolution equation can be obtained from Eqs. 101 112|) 

dv ( - .C$13 A - 
1— = H +i — - — A v . (13) 
dt \ dt ' y ' 
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where 



U ! "o Am 2 31 /2E + Vs 2 3 ) " 4) 



H 2 = ^U 12 M 2 u\ 2 + Vc 2 13 , M 2 = diag(0, Am 2 21 ). (15) 

In (|13() the matrix A has all zero elements but A13 = — A31 = 1. The last term in Eq. (|13|) can be evaluated 
as dOis/dt ~ A813/ R^th for trajectories crossing the mantle only. For the core crossing trajectories the 
derivative can be large at the border between the core and the mantle. However, because of averaging of 
oscillations driven by Am 31 we can neglect this dependence too. 

According to Eq. flU , the state v T decouples from the rest of the system and evolves independently. 
The (v e , v^) sub-system evolves according to the 2x2 Hamiltonian H 2 (v e — v^ sub-matrix in Eq. (|14H ). This 
Hamiltonian is determined by the solar oscillation parameters Am 21 , tan 2 #12 and the potential V « Vc\ 3 . 
Thus, in the propagation basis the three neutrino problem is reduced to two neutrino problem. 

Correspondingly, the evolution matrix S (the matrix of amplitudes) in the propagation basis (£ e , u T ) 
has the following form: 

(A e e A e ^ \ 

V 4 , (16) 

A TT J 

where 

A TT = cxpHcfe) , 03 « £ , (17) 

and L is the total distance traveled by neutrinos. Other amplitudes in (|16|) . A ee , j4 e „,... should be found 
by solving the two neutrino evolution equation with the Hamiltonian H 2 . 

Since the oscillations driven by Am 31 are averaged out, the dependence of the effects on the type of 
mass hierarchy (normal, inverted) appears only in the value of §13 for neutrino and antineutrino channel. 
In what follows we will present results for normal mass hierarchy. For inverted hierarchy the difference of 
results is very small and can be largely absorbed in redefinition of S13. 



B. Flavor transitions 



Let us find the probabilities of the «-> v e oscillations, P^e, and the v e «-* v e oscillations, P ee , relevant 
for our problem. The calculation proceeds in the three steps (see the transition scheme presented in fig.nj: 
(1) projection of the initial flavor state on to the propagation basis, (2) evolution in the propagation basis, 
(3) projection of the result of the evolution in the propagation basis on to the final flavor state. According 
to this picture the S— matrix in the flavor basis equals: 

S = USU\ (18) 

where U and S are given by Eqs. JT^ and (fTB|> . 
Using Eqs. fllHllflll^ we find 

+ Sl 3 £ 2 13 s 2 23 . (19) 



p — 

1 fie — 



-Sl3Cl3S23-4e 



G 



For the sub-GeV sample the oscillations driven by Am 31 are averaged out, so that there is no interference 
effect due to state v T . At the same time, according to (|19fl the amplitudes A ee and A^ e interfere. It is 
this interference which produces effect we interested in this paper. Notice that amplitudes A ee and A^ e are 
both due to the solar oscillation parameters. However, their interference appears due to presence of the 
third neutrino (non-zero si 3 ). In the limit s\ 3 = the interference disappears. In what follows we will call 
the interference of the amplitudes (with solar oscillation parameters) due to non-zero U e3 ~ S13 as the U e3 
induced interference. 

According to l|19|l . there is no interference of the amplitudes driven by the atmospheric, Am 31 , and 
solar Am^ mass splittings. This interference is averaged out for the most part of the zenith angles. If 
cos0„ > (above the horizon), neutrinos propagate in the atmosphere, where the matter effect can be 
neglected. The effect of corresponding interference terms is very small: below (0.2 - 0.3) % (see Appendix), 
though we take it into account in our numerical calculations. 

The probability (|19|l can be written explicitly as 

Pp.c = c 2 13 4 3 P 2 - 2S 13 c 2 13 S23C2 3 {cosScpR2 - sinScph) + S 2 13 c 2 13 sl 3 (2 - P 2 ), (20) 

where 

P 2 = \A^ e \ 2 = l-\A ee \ 2 , R 2 =Re(A; e A ee ), h = Im(Al e A ee ) (21) 

are the 2v probabilities in the propagation basis. 
Similarly, we get P ee : 

P ee = c| 3 (l-P 2 )+s i3 . (22) 

No induced interference appears here due to zero projection of v e on to state (see l|I2(l ). 

For antineutrinos, the probabilities P 2) R2, h should be obtained by replacement of V — > —V in the 
Hamiltonian H 2 of Eq. (fT4")l . and the sign of phase Sep should be changed: 

P 2 = P2(-V e ), R 2 = R 2 (-V e ), I 2 =h(-V e ), S la = s la (-V e ), 6 C p^-6 C p- (23) 

As a result, 

P Me = c\ 3 c\ 3 P2 - 2si 3 Ci 3 S23C23(cos<5cpP2 + sin S C p h + Si 3 C 13 S2 3 (2 - P 2 ), (24) 



P ee =4(1 ~h) +t 13 . (25) 

Let us consider the two neutrino probabilities P2, P2, I2 as well as P 2 , R2, I2 in details. We have 
calculated them numerically (see results in the figs. |3 - El using the distribution of density in the Earth 
from Ref. 0. 

Properties of the probabilities can be well understood using their expressions in medium with constant 
density: 

P 2 = sin 2 29™ sin 2 ^ , (26) 



R 2 = - sin 26>™ cos 2(9™ sin 2 ^ , (27) 
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sin26>"2sin0 m . (28) 



Here m = AH\ 2 R ear th cos 8„, is the phase of oscillations in matter, where AH\ 2 is the difference of the 
eigenvalues of the Hamiltonian in matter, R ea rth is the Earth radius and Q u is the zenith angle of neutrinos. 
In 1261 - 128[) 9^2 is the 1-2 mixing in matter determined by 

8in20K = — 2 . (29) 

cos 2 20i 2 (l - E U /E R ) 2 + sin 2 29 12 

The resonance neutrino energy equals 

Am| lC os20 12 _ n OQO n ^ T t Am 2 21 \ /2.0g/cm 3 

V7- 10- 5 eV 2 / V 57 

In the mantle, for the present best fit value Am 21 — 7.3 • 10~ 5 eV 2 and for sin 2 29i 2 — 0.8 we get 
E R = 0.10 GeV which is below the threshold of sub-GeV range. Therefore for Am^ ~ (5 — 7) • 10~ 5 
eV 2 and E v ~ (0.1 — 0.5) GeV the oscillations occur in the matter dominated regime when the potential is 
larger than the kinetic term: V > Am 2 /2E. 

For Am^/E^ < 10~ 4 eV 2 /MeV the depth of oscillations is roughly proportional to (Am 2 ) 2 . The 
oscillation length, l m , is close to the refraction length, Iq, and only weakly depends on energy: 

sin 2 20 m ^sin 2 20 12 (^^|) 2 , l m ^l = y. (31) 

With increase of Am 2 /E, the mixing parameter sin 2 20 m , and consequently, P 2 approach 1 in the 
resonance in the neutrino channel. In the antineutrino channel the mixing and P 2 increase but they are 
always below vacuum values. 

The propagation (at least in the mantle of the Earth) has a character of oscillations with quasi constant 
depth and length. Correspondingly, P 2 , R 2 and I 2 have an oscillatory behavior with cosO^. 

In Fig. [21 (upper panel) we show dependence of P2 on the zenith angle of neutrino, 0^, for different 
values of Am 21 /E. The depth of oscillation of P 2 is determined basically by sin 2 28^. P 2 monotonously 
increases with Am 21 . Notice that the first oscillation maximum is achieved at cos8„ ~ — 0.35-; — 0.4 and 
the effect is zero at cosB^ ~ —0.64. Second maximum is for the trajectories at the border between core 
and mantle: cos0„ = —0.84. For cos 8„ < —0.84 neutrinos cross both the mantle and core of the Earth. 
The interplay of the oscillations in the mantle and in the core leads to some enhancement of the transition 
probability in spite of larger density of the core. For core crossing trajectories the period of oscillation is 
smaller. 

For antineutrinos (fig- EI bottom panel) the mixing angle is suppressed. The oscillation length is smaller 
than Iq. With increase of Am 2 /E the mixing (depth of oscillations) increases whereas the oscillation length 
decreases approaching vacuum values. 

The oscillation effects in the antineutrino channel are smaller by factor 2-3. 

R 2 has similar oscillatory dependence on cosO^ (fig. |3J) with the depth of oscillations given by ~ 
sin 20J5 cos 20^ ( see Eq. I|27|)'l. In contrast to P 2 , with increase of Am 21 the real part, R 2 , first increases, 
reaches maximum at Amfi ~ 7 • 10~ 5 eV 2 (Am 21 /E u ~ 2 • 10~ 4 eV 2 /GeV ) and then decreases. The 
interference term is zero in the resonance: Am 2 2l /E v - 7 • 10~ 4 eV 2 /GeV. It changes the sign with further 
increase of Am 2 i/E v approaching vacuum value. 

In general (without rely on constant density approximation) the real part of the interference term can 
be written as 



R 2 = Re(A ee A* ) = ^/P 2 {1- P 2 )cos(^ ee - <^ e ), (32) 
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where cf> ee = arg(A ee ) and <p^ e = arg(A fie ). From Eq. (|32|) we conclude that maximal value equals 
R™ ax = |. It corresponds to P2 = 1/2 and <p ee = <j)^ e +irk, (k = integer). For the sub-GeV sample we find 
that P-2 = 1/2 is achieved at Am|i ~ 7 • 10~ 5 eV 2 , that is, for the present best fit value. 



In the constant density approximation the phase factor equals : 





1 — sin (j) m sin 29\2 



From this equation we find that in maximum of the oscillation probability (sin</> TO = 1): cos(0 ee — (p^e) = 1, 
and consequently, the interference term reaches maximum. 

According to PJ1IT7|) 



and therefore the interference probability dominates at high energies or low Am^, when 26 m < 7r/4. The 
latter corresponds to ^m% 1 /E v ~ 2 ■ 1CT 4 eV 2 /GeV for the mantle of the Earth. Thus, for Am^ ~ 7 • 1CT 5 
eV 2 , P2 and R2 are comparable. For larger Am^, the LMA probability P2 dominates, whereas for smaller 
Am^, the interference probability R2 is larger. 

The interference term R2 has opposite sign for neutrinos and antineutrinos (fig. bottom panel) due 
to change of the sign of V. This result can be easily understood using the constant density approximation 
(|27p. Indeed, the mixing angle in matter, fflfe, differs for neutrino and antineutrino. For definiteness, let us 
assume that vacuum mixing angle is below 7r/4, as is favored by the present solar neutrino data. In this case 
matter suppresses the mixing in the antineutrino channel, and enhances mixing in the neutrino channel. So, 
we have 0^{p) < 812 < 7r/4, and O^v) > #12- Furthermore, for Am 2 < 10~ 4 eV 2 (where the interference 
effect is large) and for neutrino energies relevant for the sub-GeV sample, the mixing is above resonance: 
0™ 2 { v ) > 7r /4. Therefore cos205S is positive for antineutrinos and negative for neutrinos, and since sin26*"2 
is positive in both channels the interference term has opposite sign for neutrinos and antineutrinos. 

Also behavior of the interference term R2 in the antineutrino channel with energy differs from that of R2 . 
In the antineutrino channel 20^2 ty) increases. Correspondingly, R2 reaches maximum when 29 , {^(v) = 7r/4 
{Amlx/Ey - 4 • 1CT 4 eV 2 /MeV) and then it decreases. 

The imaginary part, I2, (fig- El changes the sign with increase of the oscillation phase, and consequently, 
with cos0 2 . So, integration over the zenith angle leads to strong suppression of I2, and therefore, the CP- 
violating effects. The depth of oscillations increases according to sin20"2/2 and maximal value, 1% — 1/2, 
is achieved in the resonance, Am^i/Eu ~ 7 • 1CP 4 eV 2 /GeV. 



El 

J?2 



tan 6 



tin 



12- 



(34) 



C. Neutrino fluxes in presence of oscillations 



Let and F° be the electron and muon neutrino fluxes at the detector in the absence of oscillations. 
Then, the flux with oscillations can be written as 




where 



r{E,Q v ) 



F?(E,Q„) 



is the ratio of the original fluxes. In the sub-GeV range the ratio r depends both on the zenith angle and 
on the neutrino energy rather weakly and can be approximated by r = 2.04 — 2.06. 



Inserting the probabilities P ee and P^ e from Eqs. (1221) and (|20|) in Eq. we get expression for the 
relative change of the v e — flux: 



§ - 1 = (r4 - 1)P 2 - 
-rsi 3 Ci 3 sin 26*23 (cos (5 C p #2 - sin<5cp £2) - 

-2S 2 3 (1 - r S 2 23 ) - S 2 13 P 2 (r - 2) + 4(1 - r^ 3 )(2 - P 2 ). (36) 
Let us consider the terms of this equation in order. 

The first term on the right hand side (zero order in s 13 ) corresponds to the LMA contribution we have 
discussed in [l^. Being proportional to P2 this term increases with Am^/B up to the resonance value 
ArriQi/E = 7 • eV 2 /GeV, where P2 ~ 1. The probability is screened by the factor (rc| 3 — 1). Since 
r 2 it leads to excess of the flux for #23 < 45° and to deficiency for 6*23 > 45°. For #23 = 45° the screening 
factor equals 0.02 - 0.03. This term does not depend on S13. 

The second term in l|36[) is the effect of induced interference. It has the following properties. 

• The term depends on S13 linearly and therefore its effect may not be strongly suppressed even for 
small S13. The interference depends on the sign of S13. 

• The interference term does not have screening factor, so it can dominate for 2-3 mixing close to 
maximum. Its smallness is mainly due to smallness of S13 as well as R2 and 

• The interference term is proportional to sin 2#23 and therefore it is sensitive to the sign of #23 • 

• With increase of Am^/B the real part (similarly to R2) first increases, reaches maximum at 
A1TI21/E = 2 ■ 10~ 4 eV 2 /GeV, and then decreases and changes the sign in the resonance. The 
imaginary part increases up to the resonance value of Am|i/-E where I™ ax = 1/2. 

• For antineutrinos the interference term (as R 2 ) has the opposite sign with respect to the neutrino 
term. The amplitudes of the imaginary part have the same sign for neutrinos and antineutrinos. R 2 
reaches value 1/2 at higher Am^/E than R2 does. 

Last three terms in Eq. I|36|) are of the order s 2 3 or of higher power of §13. Practically among these terms 
only the first one can give significant contribution provided that the (2 - 3) mixing deviates from maximum. 
This term does not depend on /S.m\ 2 /E. Besides s 2 3 suppression the second term has an additional small 
factor (r—2). Its contribution does not exceed 0.1%. The third term is proportional to s 4 3 and also contains 
screening factor. 

For exactly maximal 2-3 mixing and r — 2 we get from (I36[l : 

p 

- 1 = -rsi3C 2 3 (cos(5 R 2 - sin<5 J 2 ). (37) 

That is, only the interference term gives a contribution. Since in the sub-GeV sample r = 2.04 — 2.06, no 
complete cancellation is possible. 

In what follows we will describe the deviation of the 2-3 mixing by the parameter 

^23 = \ ~ 4- (38) 
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From the 2v analysis of the atmospheric neutrino data (JIJ we get 

I-Dasl < 0.15, (90%C.L.). (39) 

Note that for consistency such a bound should be obtained from the Zv analysis which includes the LMA 
oscillations. 

III. OSCILLATION EFFECTS IN THE e-LIKE EVENTS 

In what follows we calculate dependences of the number of e— like events on the zenith angle of electron, 
e , and the electron energy. The general expression for the number of e-like events, N e , as a function of O e 
is 

N e oc V / dE u dE e d(co$Q v )dh FJE V , x 
^ J dE e 

xV(e e ,G v ,E u )K e {h,cose u ,E v )E(E e ) , (40) 

where F e is the atmospheric z^-flux at the detector given in Eq. (|35|1 (the fluxes F® and F® without 
oscillations are taken from Ref. [sEj): da/dE e are the differential cross sections taken from Ref. j^], n e 
is the normalized distribution of neutrino production points, h is the height of production, e(E e ) is the 
detection efficiency of the electron, is the "dispersion" function which describes deviation of the lepton 
zenith angle from the neutrino zenith angle ( for details see Ref. |37|). 

The integration over the neutrino zenith angle and neutrino energy leads to significant smearing of the 
O e dependence. The average angle between the neutrino and the outgoing charged lepton is about 60° in 
the sub-GeV range. Furthermore, neutrinos and antineutrinos of a given flavor are not distinguished in 
the present atmospheric neutrino experiments, so that their signals are summed in Eq. I|4U|I which leads 
typically to further weakening of the oscillation effect. 

According to Eqs. Ij^6(l and (|40|l the relative change of the e— like events, can be represented as the 
sum of three contributions: 

e e = ^-l = e L e MA + ef t + e u e «\ (41) 



where 



(rc 2 23 - 1)[(1 - 0(P 2 ) + m)} = (rD 23 + 0.5r - 1)[(1 - £)<P 2 ) + £<P 2 >] (42) 



is the contribution of oscillations driven by the solar (LMA) parameters, 

e e Tlt = cos 6 - sin 6 e\ nt , (43) 

« -r(5ia^3) «n 28 23 [(1-0 (R2) + ^(R 2 )] (44) 

e\ nt » -r(h 3 c 2 13 ) sin 29 23 [(1 - 0(h) - C<J a )] , (45) 

is the interference term, 

£e /e3 = -2(S? 3 )(l-r4) (46) 
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is the U e 3 -induced term. Here (X) = (X(E e ,Q e )} is the quantity X averaged over appropriate energy and 
zenith angle intervals of neutrino as well as final lepton (|40|) ; they are functions of the electron energy E e 
and the zenith angle O e . Here we have summed up the effect of neutrinos and antineutrinos, assuming that 
the detector does not identify the electric charge of the lepton. The parameter 

Nn 

£ = -= — — = 0.3 (47) 

describes the relative contribution of the antineutrino flux without oscillations. In estimations one can take 
(E) = 0.4 GeV as an effective neutrino energy relevant for the sub-GeV sample of events. 

Properties of different contributions H41JI reproduce basically the properties of the corresponding terms 
in the expression for the flux change l|3()[) . New features of e e are related to the strong averaging effect due 
to integration over the neutrino energy and zenith angle as well as due to summation of the neutrino and 
antineutrino signals. 

In what follows we take the experimental data, fluxes, features of detection from Ref. [l^. Recently 
Super-Kamiokande collaboration has published results of the refined analysis of the data (see e.g., Q)- In 
particular, new detector simulation, data analysis, input atmospheric neutrino fluxes and cross sections have 
been used. Unfortunately, published till now information is not enough to update our calculation. At the 
same time, we expect that impact of these changes on our results will not be significant. 



A. The LMA contribution 



Let us assume that S13 is zero or very small, so that 

e e ^e L e MA ^rD 23 {P 2 ) v „ (48) 

where {P 2 ) vB = [(1 - £){P 2 ) + £(P 2 )]. 

In Fig. we show the zenith angle dependences of the relative excess of the e-like events for different 
values of Am 2 2 . The upper panel corresponds to large deviation of the 2-3 mixing from maximum, so that 
the screening factor equals 0.33. Increase of the excess with Am 2 2 follows the dependence of (P 2 ) and 
proceeds according to increase of sin 2 29^. The zenith angle appears due to oscillations of high energy part 
of the sample. For the best fit point of the LMA solution the excess is (3 - 4)%. 

The integrated over the zenith angle excess (which corresponds to the result of the third zenith angle 
bin) can be estimated in the following way. At Am 21 ~ 7- 10~ 5 eV 2 and sin 2 26 12 ~ 0.82 the effective mixing 
parameter sin 2 26*" 2 ~ 0.35. This parameter determines the depth of oscillations of P 2 . The averaging over 
the zenith angle gives (P 2 ) ~ 0.11. (Notice that we average here over the whole interval cos0„ = — 1 -j- +1 
and the oscillation effect for the down-going neutrinos is small.) The averaging over the neutrino and 
antineutrino fluxes leads to (P 2 ) u o ~ 0.09. So, e e ss 0.09 x 0.33 = 0.03 in agreement with exact calculations. 

For large Am 2 2 the oscillations can explain the experimental results without additional rcnormalization 
of the original neutrino flux. For smaller Am 2 2 , the data points can be reproduced as a sum of the effects 
of oscillations and flux renormalization. 

For sin 2 23 = 0.45 (sin 2 28 23 = 0.99) the screening is much stronger: 0.127. Since the dependence of 
the excess on sin 2 26*23 factors out, the excess scales as the screening factor: the increase of 2-3 mixing leads 
to decrease of the excess by the 2.6 (see fig.|Hlbottom panel). This effect can be seen also in fig.Elwhere we 
show the zenith angle dependence of the ratio of events with and without oscillations. 

For sin 2 6*23 > 0.5 the oscillations produce a deficiency of the e-like events. The histograms are nearly 
mirror reflection with respect to N° sc /N® 1. According to fig.[f)]the present data disfavor values sin 2 6*23 > 
0.5 which lead to deficit of the e-like events. Thus, for Am 2 2 = 7.3 ■ 10~ 5 eV 2 we find the bound D23 < 0.1 
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(without renormalization of the original fluxes). It corresponds to a situation when all experimental points 
with error bars in fig. H3 are above the predicted curve. Let us remind that the recent cosmic ray data 
tend to decrease the original neutrino fluxes which strengthens the bound. For sin 2 #23 < 0.5 the bound on 
deviation from maximal mixing is substantially weaker: D23 < 0.4 (or sin 2 2#23 > 0.36). It corresponds to 
a situation when all experimental points with error bars in fig. are below the predicted curve. 

The dependence of the excess integrated over the zenith angle on the electron energy is shown in fig. 
The excess increases with decrease of energy according to increase of sin 29^ or P 2 (fig. [2J ■ In the very low 
energy bin, E < 0.25 GeV, the excess can reach 5 - 6% for the best fit point. The LMA contribution and 
~ 3% renormalization of the flux can give good description of the data. The excess increases from high 
energies to low energies by about 6%. Therefore measurements of the energy dependence with accuracy 
~ 2% will allow to establish existence of the LMA contribution. 

With change of the vacuum 1-2 mixing the probability P 2 changes only very weakly. 

The oscillation effect depends mainly on Amf 2 and sin 2 #23- After precise determination of the Ato 2 2 
(KamLAND will reach 10% accuracy and also SNO will contribute) one can use data on N° sc /N® in 
atmospheric neutrinos to search for deviation of 2-3 mixing from maximal. In fig. [5] we show contours 
of constant relative change of the e-like events in cos 2 #23 — Am 2 2 plane. Notice that the lines are not 
symmetric with respect to cos 2 823 = 0.5 due to deviation of r from 2. For this reason the bound from the 
side cos 2 $23 < 0.5 is stronger. We assume here that S13 w 0. Uncertainties due to unknown values of S13 
and 6c p will be discussed in sect IVC. 



B. Ue3 induced interference 

Let us assume that 1-3 mixing is non-zero but 6c p — (or tt). Now all three terms in l|41|l give 
contributions to the oscillation effect. The interference term contribution is determined by the real part R2. 
It dominates if the 2-3 mixing is close to maximal. In our estimations below we will use Am 2 2 = 7.3 • 10 -5 
eV 2 and sin 2 20 12 ~ 0.82. 

In fig. [5] (upper panel) we show the zenith angle distribution of the total oscillation effect for different 
values of S13 and s| 3 = 0.45. The LMA contribution is positive and relatively small: its value averaged 
over the zenith angle equals e^ MA = 0.8% (it is given by the line S13 = 0, see the upper panel). The U e 3 
contribution is also suppressed by the screening factor. It is negative for s 2 3 < 0.5, and being quadratic in 
S13, does not depend on the sign of 1-3 mixing. We obtain 




(Notice that there is a small matter effect on S13 which is different in neutrino and antineutrino channel.) 
The interference term is linear in S13 and positive for the negative sign of S13: 

ef* = -2.0% (|^) . (50) 

It can be estimated as follows. The depth of oscillations equals ~ sin 26*^ cos 29" 2 — 0.47. The averaging 
over cosO„ gives (R2) ~ 0.14. Due to negative effect for antineutrinos the total effect is (R2) ~ 0.09. 
Averaging over the energy (which is important here) leads to further reduction by about 30 %. Finally we 
get e™* = (i?2)^si3 ~ 2% for S13 = —0.16 in agreement with calculations. 

Summing up all the contribution we find for S13 = —0.16: e* ot = 2.5% in agreement with result 
in fig. El (upper panel). For S13 = +0.16, the interference term changes the sign and we get: e' ot = 
— 1.7%. Apparently the curves are not symmetric with respect to N° sc /N® = 1 due to the LMA- and f7 e 3- 
contributions which do not change the sign with si 3 . When |si 3 | decreases, both e™* and e^ e3 decrease in 
absolute value. 
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In fig. [5] (lower panel) we show the zenith angle dependences for larger deviation of 2-3 mixing from 
maximal value, = 0.35. Now screening is weaker and the LMA oscillations give main contribution 
€ lma = 2 9% This leads t0 the smft f all th e histograms to N° sc /N® > 1. 

Now also the U e % contribution is larger being comparable with the interference contribution: 

2 



e Ue3 « -1 



In contrast, the interference term has no screening factor and its absolute value even slightly decreases in 
comparison with the previous case due to decrease of sin 2 2#23 : 



• 9% (S)- (52) 



This leads to more complicated dependence of the excess on s\ 3 . We find that maximum of total excess is 
realized for S13 = —0.16: e* ot = 3.6%. For si 3 = —0.08 we get a very similar value: e* ot = 3.55%. 

Maximal effect of interference can be estimated in the following way. As we have discussed, R™ ax = 1/2 
(which can be achieved at Am 2 2 = 7.3 • 10~ 5 eV 2 ). The averaging over the zenith angle gives (R2) = 
R™ ax /4 = 0.125. The antineutrino contribution is negative, so that ef* < 2s 13 (l - £)0.125 = 0.18si 3 < 3%. 
Averaging over the energy leads to an additional suppression. 

Using curves which correspond to different sign of S13, it is easy to disentangle contribution from the 
interference term. Obviously, 

C t M = lK° t (s 13 )-ei° t (-s 13 )], (53) 
and for the two other contributions we get: 

e LMA = e tot {si3 = 0)j e Ue3 = 1 [ e ft M + e ^(_ Sl3 )] _ e lot {si3 = Q). ( 54 ) 

In Fig. we show dependence of the ratio of number of events integrated over the zenith angle, 
N° sc /N®, on the energy. According to our analytical consideration, with decrease of energy the LMA 
contribution (green histogram) increases fast, the C/ e 3-contribution is unchanged and the interference term 
first, increases but then below E ~ 0.4 GeV starts to decrease. Using relations (|53l - 154|) we find from the 
Fig. H@ (upper panel) for s vi = -0.16 that in the bins E = (0.4 - 0.65), (0.25 - 0.40), (0.10 - 0.25) GeV: 
€ lma = L2%i 2,7%, 5.9% respectively, ef e3 = -1.1% = const, whereas ef* = 1-7%, 2.7%, 2.4%. At 
high energies the interference term gives main contribution. (Notice, however, that for E > 0.6 GeV our 
approximation may not be precise). 



C. CP-violation effects 



Let us consider effects of the CP-violating phase 6c p- Notice that if we substitute I2 by its vacuum 
value (eq. (|20|) with 6 m — > 6), the interference term, as is expected, becomes equal e* n * = 1/2 AP, where 
AP is the neutrino-antineutrino CP- asymmetry. 

According to fig.^ I2 alternates the sign with change of the zenith angle. However, there is no averaging 
of the effect to zero for two reasons: 

1. For the mantle trajectories up to 1.5 - 2 periods of oscillations are obtained. In particular, for 2 • 10~ 4 
eV 2 /GeV (fig. 0J| which corresponds to the best value of Am 2 2 and E = 0.4 GeV there are 1.5 periods. 

2. Due to change of the density for trajectories with different e the curves are not symmetric with 
respect to I 2 = 0. 
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For antineutrinos the probability I2 has smaller amplitude and oscillation length, furthermore, the 
curves are nearly symmetric with respect to I2 = 0. As a results, integration over the zenith angles 
leads to strong suppression of the averaged value (12). This, as well as difference of the original neutrino 
and antineutrino fluxes result in existence of the CP-odd effects, even in the sample where neutrino and 
antineutrino signals are summed up. 

In fig. El we show the zenith angle dependence of the oscillation effect for different values of Sep- The 
analytical expression of this dependence on Sep is given in Eq. (|37[) . Simple estimation of the effect can 
be obtained as follows. Using results of fig. 0] we obtain after averaging over the zenith angle and energies: 
(I2) = 0.040 and (I2) ~ 0.018. Then for S13 = —0.16 and s| 3 = 0.35, the contribution of imaginary part to 
the interference term equals: e™* ~ 1.1%. In the previous section we have found that for the same set of 
parameters e^* w 2.0%. So, the interference effect can be written as 

4 nt « (2.0cos<5 C p-l.lsin£op)% (4s = 0.35). (55) 

The relative values of numerical coefficients depend on the 2v probabilities R2 and 12- For 5c p — — 7r/4 
we find 4 n * « 2.2% (no change with respect to 8 C p = 0). For S C p = tt/2: 4™* = 4"* ~ Maximal 
effect of I2 is for Sep = —tt/2: e* nt = e™* s» 1.1%, that is, the excess decreases by 2.8% in comparison with 
5 C p = case (see fig. ITT1 upper panel). From fig. Eland relations (|531 - we find e^ MA + e^ e3 w 2.2%, 
and consequently, 

4 ot = 2.2% + e™*. (56) 

The latter formula reproduces well the results shown in fig. ^2 

Let us compare different contributions to the excess of the e-like events. As we have established the 
largest contribution can be obtained from the LMA term. Maximum of e^ MA ~ 6% is achieved for the largest 
possible Amj 2 , minimal energy and largest deviation of the 2-3 mixing from maximum. The contribution 
does not depend on S13. 

The interference term gives maximal contribution, e™* ~ 3%, for Ato 2 2 = 7 • 10" 5 eV 2 and maximal 
possible value S13. It depends very weakly on the deviation -D23- 

The U e 3 maximal contribution, e^ e3 ~ 1 — 2%, is realized for the largest possible values of S13 and Z?23- 
It does not depend on Am 2 2 in the first approximation. According to l|36|) this term has an opposite sign 
with respect to the LMA term and therefore partial cancellation with the LMA contribution always occurs. 

These results allow one to understand that the contributions of non-zero S13 (interference term and 
U e z) can not further enhance the excess produced by the LMA term. Indeed, for large Am 2 2 ~ 2 • 10~ 4 eV 2 
where e^ MA is maximal, the interference term contribution is already small, and moreover, the U e % term is 
negative compensating substantially the positive e*™* contribution. As a result the LMA contribution can 
be enhanced by about 1% at most. For Am 2 2 ~ 7 • 10~ 5 eV 2 where the interference term is maximal, the 
LMA contribution is smaller and again partial cancellation with U e ^ contribution occurs. 

Notice that the cancellation of the LMA induced excess can be stronger than the enhancement since 
e* n * and e^ e3 can have both the same negative sign with respect to e^ MA . 



IV. MEASURING D 23 AND 5 C p 

In terms of the deviation parameter l|38ll the excess of the e-like events can be written as 

e e = r «P 2 ) - 2S 2 3 ) D 23 + ef* + (r/2 - 1) «P 2 ) + 2S 2 3 ) . (57) 

The interference term depends on the deviation very weakly: e*™* cx sin 26*23 = V 1 - 4£) 23 « 1 " 2D h- 
Since variations of D23 in the presently allowed range (|3*^| change e™* by less than 5% and in the first 
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approximation this dependence can be neglected. Also the last term in (|57|> does not exceed 0.5 %. Therefore 
with a good approximation e e is a linear function of Z?23- We can find coefficients of this function using 

% El 

For Am 2 2 ~ 7 • 10~ 5 eV 2 and zero 1-3 mixing we get 

e° e » (20.5D 23 - 0.3)%. (58) 

For | s 13 1 = 0.16 maximal (corresponds to S13 = —0.16) and minimal (S13 = 0.16) values of e e can be 
approximated as 

e j«« = (i2.5Da3 + 1.8)%, C m = (12.8^23-2.1)%. (59) 

Since for the present upper limit on D23, e™ m (0.15) < 0, any upper experimental bound, ef° p , will not 
improve the limit for Z? 23 l|39|) . 

If the lower bound on e e is established, one can put the lower bound on -D 23 . Using expression for 
e rnax we find tllat the lower bounds £ ezp > 2%j and £ exp > 3% will give ^ > q.02 and L> 23 > 0.10 

correspondingly. 

The bound on £) 23 can be improved if future experiments put stronger bound on S13. For | S13 j = 0.08 
maximal and minimal values of e e equal 

e max w ( 17£>23 + .84)%, e™ m « (18.5L> 23 - 1.2)%. (60) 

From the expression for e™ m we find that the present bound on Z? 23 will be improved provided that the 
upper bound on the excess is better than 1.6 %. If e e e xp < 1%, we get D 23 < 0.1 etc.. 

Using e™ ax we obtain that the lower experimental bound e e e xp > 2% will lead to the lower bound 
D23 > 0.07. 

Suppose very strong bound on S13 is obtained and the lower bound e exp > 0.02 will be established. 
Then according to fig.[S]the interval Am 2 2 = (7.3 ±0.7) • 10" 5 eV 2 (10% error) will lead to the lower bound 
on deviation from maximal mixing D23 > 0.1. 



Can the CP-violation phase be measured? As follows from the fig. ^2 the phase Sep does not produce 
any particular zenith angle dependence and the energy dependence (not showed). The same effect can be 
achieved by changing other parameters. 

Let us consider the most favorable case: Am 2 2 ~ 7.3 ■ 10~ 5 eV 2 and \sis\ ~ 0.16. The total relative 
change of number of the e-like events can be written as 

e e = (l.ScosScp - 0.8sin£ CP ) + 13.6£> 23 %. (61) 

Depending on Sep the interference term changes in the limits —1.92 +1.92. The LMA contribution is 
restricted by \e^ MA \ w 13.6|£> 23 | < 2.04%. So, the predictions can be in the interval -3.96 +3.96. 

On the other hand, for 5cp = 0, |e™*| = 1.8%. Therefore for zero CP-violating phase the excess 
(deficiency) can be in the intervals: (0.2 3.8)% and (—3.8 -. — 0.2)% which covers practically whole the 
interval predicted for non-zero CP-violating phase. So, to get any information about 5cp one needs to 
improve the bound on the deviation D23 from independent measurements. From fig. it follows that 

A £e = -2%. A(si Q n2 / 23) . (62) 

E.g., 1% effect of Sep can be produced also by 0.05 change of sin 2 623: from 0.35 to 0.40. Variations by 3% 
would require the increase of sin 2 Q23 from 0.35 to 0.5. 
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Apparently, the effect similar to that of 5c p can be produced also by small variations of s%s, or Am^. 
Let us consider this degeneracy of parameters. 

Using fig. (upper panel) we find that change of the excess by 1% can be achieved by (35 - 40)% change 
of Amf 2 near the best fit point, e.g., from 5.2 to 8.7 • 10~ 5 eV 2 : 

Ae e = 0.44% • (63) 
10- 5 eV 2 V ' 

(for sin?^ = 0.35). The effect is smaller if the 2-3 mixing is closer to maximum. Ae e = 1% would require 
rather large change of Am 2 2 : from 5 to 15 of 10 -5 eV 2 . 

Further operation of the KamLAND and SNO will allow to determine Am 2 2 with 10% ambiguity, which 
is transferred in 0.3% ambiguity in e e . 

The degeneracy of the 5c p and S13 is more complicated and it depends on specific value and sign of 
S13 as well as sin 2 #23- In particular, for sin 2 $23 = 0.35 and S13 = —0.16 the dependence of e e on S13 is 
rather weak (fig. EJ): A reduction of the excess by 1% requires increase of S13 from - 0.16 to + 0.05. For 
sin 2 6*23 = 0.45, Ae e = 1% can be compensated by changes of S13 in the interval —0.22 -= — 0.10. That is, 
moderate accuracy of measurements of S13 could be enough to determine 5c p- Notice, however, that with 
decrease of S13 the effect of 5c p decreases. 

Dependence of the oscillation effect on 1-2 mixing is very weak: variations of tan 2 O12 in the interval 
from 0.30 to 0.52 produce a change Ae e = 0.3%. Expected improvements of determination of tan 2 d\i will 
further reduce this ambiguity. 

The main problem is the identification and measurement (or restriction) of the oscillation effect in view 
of large present uncertainties in the original neutrino flux (15 - 20 %). In principle, if high enough statistics 
will be achieved the oscillation effect can be distinguished from the renormalization by its zenith angle and 
energy dependences. At the same time, further improvements in the calculations of the neutrino fluxes are 
extremely important. Also separate measurements of the neutrino and antineutrino signals will help. 



V. CONCLUSION 



1 . After confirmation of the LMA MSW solution of the solar neutrino problem it is clear that the effect of 
v e - oscillations should appear in the atmospheric neutrinos at some level even for zero value of S13. 

For the allowed values of the oscillation parameters, in particular, sin 2 26*23 > 0.91, the LMA oscillations 
can produce the integrated effect (excess or deficit) up to (5 - 6) % in the sub-GeV sample. The effect 
increases with decrease of energy and in the low energy part of the sample it can be as large as 8%. The 
zenith angle dependence of the effect is rather weak with maximum achieved in the upward-going bins. 

The LMA effect is strongly suppressed for exactly maximal 2-3 mixing, so searches for the oscillation 
effect in the sub-GeV sample can be used to measure the deviation of 2 - 3 mixing from maximum. Here, 
however, an ambiguity appears due to unknown values of S13 and 5c p. The ambiguity can be reduced if 
stronger bound on S13 will be established from independent measurements. 

2. The present experimental accuracy is comparable with the maximal expected effect. Notice that without 
additional renormalization of the original neutrino fluxes, the data show some excess of the e-like events 
which can be explained (at least partially) by the LMA-oscillations. In fact, the data (including weak 
zenith angle and energy dependences of the excess can be perfectly reproduced by the LMA-contribution 
corresponding to the best fit values of parameters and partial (3- 5%) renormalization of spectrum. 

The excess of e-like events in the sub-GeV sample and the absence of the excess in the multi-GeV range 
(as it is indicated by the present data) testify for the deviation of the 2-3 mixing from maximum. 
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3. Non-zero 1-3 mixing gives an additional contribution to the oscillation effect. For the sub-GeV sample, 
it leads to interference of the the two neutrino amplitudes driven by solar oscillation parameters (induced 
interference). The interference term is linear in S13 and does not contain the screening factor. It dominates 
if 2-3 mixing is close to maximal. 

The interference term can reach 2 - 3%. The maximal value (real part) corresponds to Am 2 2 w 7 • 1CP 5 
eV 2 for the sub-GeV sample. 

The 1-3 mixing leads also to contribution proportional to s 2 3 which does not exceed ~ 1%. 

4. The interference term depends on the CP-violating phase 5c p- Variations of 5c p can change the 
oscillation effect by |Ae| = 3%. 

5. The relative effects of 5c p is enhanced for the sample induced by neutrinos or antincutrinos, that is, 
when the sign of the electric charge of electron is identified. 

6. Various contributions to the oscillation effect have slightly different zenith angle and also energy depen- 
dences. This, in principle, can be used to separate them. In particular, the LMA-contribution increases 
with decrease of the energy, the interference term first increases and then starts to decrease. 

The zenith angle dependence is very weak for low energy bins. 

7. There is strong "degeneracy" of parameters once total excess is measured only. The same integral 
oscillation effect can be produced for different values of sin 2 823, Am 2 2 , S13 and Sep- 

8. In principle, future high statistics studies of the atmospheric neutrinos will allow to measure the neutrino 
oscillation parameters. For this the accuracy of measurements of the oscillation effect should be about 1% 
or better. Also a way should be found to distinguish the oscillation effect from the effect of the neutrino 
flux normalization. The problem of degeneracy of parameters should be resolved. There is a good chance to 
measure Am 2 2 with high enough accuracy, so that the corresponding uncertainty will be eliminated. It will 
be very difficult to resolve ambiguity related to of sin 2 #23, S13 and Sep- If future (e.g., reactor) experiments 
put stronger bound on S13, the ambiguity related to s 13 and 5c p can be substantially reduced. This will 
allow to use the atmospheric data to restrict a deviation of the 2-3 mixing from the maximal one. 

9. With present knowledge of the oscillation parameters, one can expect the effect of v e oscillations at the 
level of existing experimental error bars and uncertainties in the normalization of fluxes. The effect of the 
LMA oscillations should be taken into account in the analysis of the atmospheric neutrino data. 
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Appendix 

Let us evaluate the effect of the interference between the solar and the atmospheric frequencies we have 
neglected in our consideration. This interference gives additional terms to the probabilities P^ e (|19f) and 

Pee G2|: 



AP F = -2s 2 3 C 2 3S23-ftr ee + 2si 3 C 2 3 S 2 3C23^e: 



(64) 



18 



(65) 



where 



K R 



Re 



A* ee A TT 



Kf,,, = Re 



'■$C1 



(66) 



Notice that K^ e depends on the CP-violating phase. Inserting (|64|l and (|65|) in l|35|) we get the corresponding 
corrections to the relative change of number of the e-like events: 



Ae P = 2t 



13 c 13 



(1 



""23 



){K ee ) + 2rsi 3 cl 3 s 2 3C 2 3{K ! 



fie/ 



(67) 



Here (...) denotes the averaging over the energy and the zenith angle. 
Let us evaluate two terms in this equation in order. 

1) . The first term is proportional to two small factors: sf 3 I?23 < 0.015. For trajectories with cos 9 > 
taking A ee k 1 we estimate: (K ee ) < (cos $3). Since for typical energy 0.4 GeV, the oscillation length in 
vacuum is I13 ~ 500 km, and the phase equals $3 ~ 2n. Therefore averaging over the zenith angle and the 
energy lead to strong suppression: (K ee ) ~ 0.2. As a result, the whole term is smaller than 0.3%. 

2) . In the second term of (|67|l {K^ e ) can be estimated in the following way. In vacuum: 



Afie — 



-S12C12 (l 



-i<t>2 



Amf 2 L 
2E 



For trajectories with cos 8 > the phase driven by the solar mass split is small: <p2 < 0.2, so that 

,2nd 1 



-Sl2Ci 2 (4>2 sin(03 + Sep)) 



-Sl2Cl2( 



l\2 COS 6 



+ S C p)), 



(68) 



(69) 



where d ~ 20 km is the depth of the atmosphere. For E — 0.4 GeV the oscillation length equals l\ 2 — 1.4- 10 4 
km. The averaging over the zenith angle gives (1/ cosO) = 3.2. Then taking sin(<^3 + Sep) — 1 we obtain 

from dnni 



\(K 



lie/ 



< 3- 10- 2 si 2 c 12 , 



(70) 



and consequently, for sxa < 0.16 the contribution of the second term to Ae e < 0.5%. Averaging over the 
energy leads to further suppression of this contribution. 
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FIG. 1: The neutrino transition scheme. Initial and final flavor states are shown in circles. In boxes we show the 
states of the propagation basis. Lines connect states between which the transitions can occur. The lines with arrows 
indicate transitions and projections relevant for oscillation channels of interest. 
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FIG. 2: Zenith angle dependences of the transition probabilities for neutrinos, P2, (upper panel), and for antineu- 
trinos (lower panel), P2, for different values of Aml 2 /E and sin 2 2#i2 = 0.82. 
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FIG. 3: Zenith angle dependences of the real part of the interference probability for neutrinos, R2 (upper panel), 
and for antineutrinos, R2 (lower panel), for different values of Am? 2 /£ and sin 2 2#i2 = 0.82. 
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FIG. 4: Zenith angle dependences of the imaginary part of the interference probability for neutrinos, I2 (upper 
panel), and for antineutrinos, I2 (lower panel), for different values of Am? 2 /£ and sin 2 2#i2 = 0.82. 
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FIG. 5: Zenith angle distributions of the e-like sub-GeV events for sin #13 = 0. The ratio of the numbers of events 
with and without i^-oscillations, N° ac /N®, for different values of Am 12 . From the upper to lower histograms the 
corresponding value of Amf 2 decreases. We take sin 2 2#i2 = 0.82 and sin 2 623 — 0.35 (upper panel), sin 2 623 = 0.45 
(lower panel). Shown are also the Super-Kamiokande experimental points from Ref. |18j ]. 
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FIG. 6: Zenith angle distributions of the e-like sub-GeV events for sin #13 = 0. The ratio of numbers of events 
with and without i/ e -oscillations, N° sc /N^, for different values of sin 2 6*23. We take sin 2 26\2 = 0.82 and Am 2 2 = 
7.3 -10" 5 eV 2 . From the upper to lower histograms the corresponding value of sin 2 623 increases. Shown are also the 
Super-Kamiokande experimental points from Ref. flftj |. 



O d> 



O" — 

0) 

O Q) 



1.15 



Am% 1= 7.3E-05 eV 



Am! 1= 6.0E-05 eV 2 



Ami 1= 4.7E-05 eV 



Am 2 21 = 20.0E-05 eV 2 
Am 2 21 = 14.5E-05 eV 2 
Am 2 21 = 10.0E-05 eV 2 



1.1 



1.05 - 




20 



FIG. 7: The excess of the e- like events integrated over the zenith angles as a function of energy. Dependence of 
the ratio of numbers of the e-like sub-GeV events with and without i/ e -oscillations, N° BC /N®, on the visible energy 
for different values of Am^. From the upper to lower histograms the corresponding value of Am 2 2 decreases. For 
other parameters we take: sin 2 2#i2 = 0.82, sin 2 623 = 0.35 and sin #13 = 0. Shown are also the Super-Kamiokande 
experimental points from Ref. |18f . 
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FIG. 9: Zenith angle distributions of the sub-GeV e-like events for non-zero 1-3 mixing. The ratio of the numbers of 
events with and without v e oscillations for different values of sin #13. We take sin 2 2(?i2 = 0.82, Am 2 2 = 7.3 ■ 10 -5 eV 2 
and sin 2 #23 = 0.45 (upper panel), and sin 2 #23 = 0.35 (lower panel). The excess increases with decrease of sin #13. 
Shown are also the Super-Kamiokande experimental points from Ref. [Tftj |. 
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FIG. 10: The excess of the e- like events integrated over the zenith angles as a function of energy. Dependence 
of the ratio of number of e-like events with oscillations and without oscillations on the visible energy for different 
values of sin #13. For other parameters we take: sin 2#23 = 0.91 and Am 12 = 7.3 • 10 eV . The excess increases 
with decrease of sin #13. Shown are also the Super-Kamiokande experimental points from Ref. |Tsj |. 
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FIG. 11: Zenith angle distributions of the sub-GeV e-like events for non-zero 1-3 mixing. The ratio of the numbers 
of events with and without v £ oscillations for different values of the CP-violating phase Sep- We take sin 2 623 = 0.35, 
sin 2 20i2 = 0.82, sin 613 = 0.16, Am 2 2 = 7.3 • 10 -5 eV 2 . In the upper (lower) panel, the excess decreases (increases) 
with increase of Sep- Shown are also the Super-Kamiokande experimental points from Ref. [18|. 



